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ABSTRACT 


Preliminary  absorptivity  data  were  taken  with  the  spectrophone 
method  and  a  1.75  watt  C02  laser.  Test  gases  included  pure  C02;  a  mixture  of 
Co2,  N  and  He;  and  a  sample  of  contaminated  He.  The  measured  absorptivity 
of  C02  at  10.6  microns  was  6.8  x  10  ^cm  This  value  agrees  closely  with  the 
absorptivity  measured  by  conventional  spectroscopy  methods.  The  absorptivity 
of  the  contaminated  He  was  on  the  order  of  5.4  x  10  cm  ^ . 

I.i  the  spectrophone  method;  a  laser  illuminates  a  closed  sample 
chamber  with  end  windows  appropriate  to  the  type  of  source  laser.  Any  energy 
absorbed  by  the  gas  causes  a  temperature  and  pressure  rise  in  the  cell.  Pressure 
rises  as  small  as  0.02  dyne  cm'2  may  be  observed  with  a  pressure-measuring 
instrument  such  as  a  sealed  capacitor  microphone.  In  the  experiment  described 
above,  the  sample  tube  was  10  cm  long.  The  absorption  pressure  signal  from 
pure  C02  was  630  dynes  cm"2.  However,  there  is  a  spurious  long-term  pressure 
rise  of  about  1  dyne  cm’2sec’1.  These  results  indicate  that  absorptivities  as 
weak  as  10'8cm_1  should  be  measurable  when  the  spurious  pressure  rise  is 
eliminated. 


Hie  spectrophone  method  with  a  ruby  laser  was  used  to  obtain  a 
spectrum  of  the  water  vapor  absorption  line  tit  6943. 8A.  The  sample  was  air 
at  one  atmosphere  saturated  with  water  vapor  at  room  temperature.  The  peak 
absorptivity  recorded  was  3  x  10  cm  .  The  linewidth  was  less  than  0.2A. 
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INTRODUCTION 

The  goal  of  this  project  is  to  measure  very  small  absorptivities 
of  light  in  gases.  The  method  capitalises  on  the  high  monochromatic  power 
available  from  lasers.  When  a  laser  beam  passes  through  ar.  absorbing  gas, 
the  absorbed  energy  causes  a  temperature  and  pressure  rise.  Pressure  rises 
as  small  as  0.02^  Hg  can  be  measured  without  difficulty  in  low  bandwidths. 

The  gas  absorptivity  is  proportional  to  the  pressure  rise.  Ttie  measure nent 
is  not  disturbed  by  a  small  amount  of  scattering  in  the  gas. 

An  earlier  report*  described  experiments  to  measure  gas  absorp¬ 
tivity  by  what  is  called  the  acoustic  method.  The  present  report  describes 
initial  experiments  using  the  spectrophone  method. 

The  spectrophone  method  is  illustrated  in  Figure  1.  The  gas 
sample  is  enclosed  in  a  cylindrical  sample  cell  with  end  windows  whose  trans¬ 
mission  characteristic  is  appropriate  to  the  source  laser.  The  cell  radius 
fits  the  beam  radius  more  or  less  closely,  but  not  so  closely  that  laser  light 
incident  on  the  tube  sidewalls  will  cause  strong  desorption.  A  short  duct 
connects  the  sample  chamber  to  the  pressure  transducer.  We  are  using  the 
MRS  Baratron,  which  is  a  sealed  capacitor  microphone.  For  a  pulsed  laser, 
the  measuring  bandwidth  is  limited  by  the  pressure  decay  time  constant  ar 
the  heat  is  dissipated  in  the  tube  walls.  For  a  CW  laser,  a  square-wave 
chopping  frequency  may  be  chosen  which  matches  the  time  constant  for  attaining 
the  steady-state  pressure.  The  energy  or  power,  and  the  laser  wavelength, 
are  recorded  along  with  the  pressure  signal. 

The  spectrophone  method  has  been  applied  successfully  with  both 
the  CO-  laser  at  10. 6n,  and  the  ruby  laser.  Aosorptivities  as  weak  as 
10  cm-1  have  been  measured.  The  available  sensitivity  appears  to  be  much 
better  even  than  this. 
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EXPERIMENTS  WITH  THE  C0„  LASER  AND  SPECTROPHONE  METHOD 

Apparatus  was  assembled  from  available  parts  to  test  the  feasibility 
of  the  C02  laser-excited  spectrophone .  The  laser  used  was  a  2-meter  long, 
nominal  4-watt  C02  laser  operating  at  10.6  microns.  The  spot  size  was  about 
1/2  cm  in  radius.  No  mode  control  was  used.  The  sample  cell  was  10  cm  long 
by  2  cm  in  radius,  with  NaC l  end  windows  perpendicular  to  the  laser  beam.  The 
cell  v/as  filled  by  evacuating  it  to  100n  Hg,  and  then  allowing  the  sample  to 
fill  it  to  atmospheric  pressure.  The  pressure  transducer  used  was  the  MKS 
Baratron,  with  a  maximum  linear  range  of  ±  1  mm  Hg.  An  oscilloscope  and 
camera  recorded  the  pressure  variations  about  the  nominal  atmospheric  pressure. 
The  laser-beam  chopper  had  a  50-percent  duty  cycle,  and  its  frequency  was 
adjustable  to  as  low  as  0.125  Hz.  Power  was  measured  with  a  crude  water 

calorimeter. 


EXPERIMENTAL  DATA 

Some  of  the  experimental  results  are  shown  in  Figure  2.  The 
gases  tested  were  pure  C02;  laser  mix,  consisting  of  5  percent  C02,  15  percent 
N  ,  and  80  percent  He;  and  He  with  an  unknown  contaminant.  The  absorptivity 
of  the  pure  CO,  was  so  strong  that  the  laser  power  input  had  to  be  reduced  to 
1.75  watts  in  order  to  observe  the  complete  response  within  the  linear  range 
of  the  Baratron.  This  reduced  power  was  used  for  all  four  photographs  in  the 

figure. 


The  step  response  of  the  pure  C02  exhibits  an  interesting  break 
point  at  0.85  mm  Hg  pressure.  The  same  break  point  is  observed  three  times 
in  the  square-wave  response,  Figure  2b.  We  demonstrated  that  this  break  point 
is  not  of  instrumental  origin.  It  is  observed  when  the  response  is  started 
in  different  portions  of  the  Baratron  linear  range.  Also,  it  is  not  observed 
on  the  cooling  cycle  of  the  pure  C02  square-wave  response.  The  break  point 
is  due  to  the  sudden  onset  of  convection  which  tends  to  cool  the  gas  more 
rapidly  than  heat  conduction  alone.  The  Rayleigh  number  for  C02  between 
parallel  plates  separated  by  2  cm  is  0.1  percent.  This  means  a  pressure 
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(c)  Laser  Mix  Square  Wave  Response  (d)  Contaminated  Helium  Square  Wave  Response 

57,  C0o.  157,  N_,  807,  He 


PERKIN  ELMER 


Report  No.  8670 


differential  of  0.76  mm  Hg  must  exist  between  the  lower  and  upper  plates  before 
convection  can  occur.  A  similar  situation  occurs  in  the  spectrophone  cell.  The 
beam  passes  through  the  center  of  the  cell.  A  temperature  gradient  exicts 
between  the  heated  gas  in  the  beam  and  the  cool  gas  near  the  top  of  the  cell. 
Since  the  Rayleigh  number  goes  inversely  as  the  cube  of  the  plate  separation, 
we  can  eliminate  this  convection  effect  in  the  next  apparatus  by  making  the 
sample  cell  fit  the  beam  more  closely. 

Figure  3  is  a  plot  of  the  calculated  step  response  from  Appendix  A. 
If  we  assume  that  convection  is  not  an  important  source  of  heat  loss  compared 
with  conduction,  we  may  use  Equation  (A-2)  to  obtain  the  absorptivity.  If  K 

for  CCL  is  1.4  x  10-4  watt/cm°C,  we  obtain  an  absorptivity  a  of  6.84  x  10  4cm  ^ . 

1  3 

This  value  agrees  closely  with  the  absorptivity  measured  by  Barker  and  Adel 

in  a  700  cm  cell  using  conventional  spectroscopy  methods. 

The  laser  mix  square-wave  response,  Figure  2c,  shows  a  much  smaller 
time  constant  because  of  the  high  conductivity  and  diffusivity  of  helium.  The 
response  satisfies  the  equation 

p  =  62. 5p  Hg  [l  -  exp  (-t/0.447  sec)] 

Figure  2d  shows  the  square-wave  response  of  contaminated  helium. 

The  response  satisfies  the  equation 

p  =  1.28p  Hg  [1-exp  (-t/0.327  sec)]  +  1.07p  Hg  t/sec. 

The  response  appears  to  be  a  combination  of  the  exponential  heating  response 
of  the  helium,  with  a  short  time  constant,  plus  a  heating  effect  with  a  very 
long  time  constant.  In  a  separate  experiment,  th.'.s  longer  time  constant  was 
observed  to  be  greater  than  several  minutes.  It  is  not  heating  and  expansion 
of  the  end  windows  or  the  cell  walls,  as  the  response  would  then  be  several 
orders  of  magnitude  below  the  observed  response.  We  are  investigating  the 
source  of  this  effect.  Assuming  the  conductivity  of  the  helium  is  14.2  x 
10"4  watt/cm°C,  the  absorptivity  of  the  contaminant  is  5.37  x  10“^cm 
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Figure  3.  Absorption  Cell  Pressure  -  Step  Response 
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SUCCESSFUL  OBSERVATION  OF  A  WEAK  WATER  VAPOR  ABSORPTION  LINE 

The  first  absorption  spectrum  successfully  taken  via  the  spectro- 

phone  method  is  shown  in  Figure  4.  The  sample  was  air  at  one  atmosphere 

pressure  saturated  with  water  vapor  at  room  temperature.  The  absorption  line 

0  -6  -1 
shown  is  centered  at  6943. 8A.  The  peak  absorptivity  is  about  3  x  10  cm  . 

This  absorptivity  was  measured  in  a  sample  cell  only  45  1/2  cm  long.  It  repre¬ 
sents  true  absorption,  not  absorption  plus  scattering.  The  absorption  line  is 
probably  somewhat  narrower  than  it  appears  in  the  figure.  It  appears  to  have 

O 

a  0.2A  width  in  th*>  figure  because  the  ruby  laser  was  operated  with  multiple 
modes,  and  because  of  pressure  broadening.  The  laser  pulse  wavelength  was 
scanned  over  a  range  of  about  1-1/2  wave  numbers  by  temperature  tuning  from 
16.3°C  to  37°C.  (In  this  report,  'wave  number"  is  defined  as  l/\.) 

EXPERIMENTAL  EQUIPMENT 

The  ruby  laser  is  shown  in  Figure  5»  (The  dye  cell  was  not  present 
in  the  cavity  for  this  experiment.)  The  ruby  rod  itself  measures  1/4  inch  in 
diameter  by  3  inches  in  length.  It  is  enclosed  in  a  coolant  Jacket  through 
which  a  mixture  of  methanol  and  water  or  ethylene  glycol  and  water  is  circulated 
to  maintain  constant  temperature  and  to  provide  temperature  tuning.  The  helical 
xenon  lamp  may  be  flashed  every  minute  and  20  seconds .  Pulses  of  300  to  400 
milli joules  may  be  obtained  with  about  10-  to  15-percent  repeatability  without 
Q-switching.  A  calorimeter  monitors  the  energy  of  every  pulse.  In  this  experi¬ 
ment  the  laser  was  fired  by  a  delay  system  which  triggered  the  recording  oscillo¬ 
scope  sweep  5  milliseconds  before  the  flashlamp.  Laser  pulses  of  1-millisecond 
duration  time  occurred  about  1  millisecond  after  the  flashlamp  was  triggered. 

THE  SAMPLE  CHAMBER 

The  sanqplc  chamber  is  shown  in  Figure  6.  The  laser  pulse  passes 
through  the  horizontal  glass  tube  shown  in  the  center  of  the  tank.  The  tube 
diameter  is  4-l/2mm.  The  laser  beam  is  limited  to  a  diameter  of  3a*n  by  an  iris 
diaphragm.  Careful  alignment  prevented  the  main  beam  from  illuminating  the  sides 
of  the  tube.  The  absence  of  the  strong  pressure  signal  with  a  1-1/2-second  decay 
constant,  which  was  previously  observed  and  reported,  is  attributed  to  thla 
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Figure  4 .  HjO  Vapor  Absorption  Line  at  6943. SoA 
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Figure  5.  Temperature-Controlled  Ruby  Laser 
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alignment  and  the  3n«n-diameter  .top.  Apparently  the  desorption  of  adsorbed 
gas  films  from  the  end  windows  of  the  sample  chamber  is  less  serious  a  prob¬ 
lem  than  desorption  from  the  chamber  sidewalls. 

The  sample  tube  is  housed  in  the  acoustic  pressure  tank  left  over 
from  experiments  evaluating  the  acoustic  method.  When  this  tank  is  closed  it 
provides  approximately  30  db  shielding  from  acoustic  room  noise.  It  also  lm- 
proves  the  thermal  stability  of  the  sample  chamber. 

The  Baratron  instrument  measures  the  differential  pressure  change 
between  the  sample  chamber  and  the  surrounding  acoustic  tank.  The  instrument 
is  attached  by  ducts  in  the  back  of  the  chamber  (not  shown  in  the  figure). 

Heat  energy  absorbed  in  the  sample  is  transmitted  to  the  sample 
tube  walls  with  a  time  constant  of  30  to  50  milliseconds.  During  this  time 
the  pressure  difference  is  observed.  In  order  to  prevent  a  pressure  buildup 
between  the  sample  chamber  and  the  surrounding  acoustic  chamber,  the  sample 
chamber  is  provided  with  a  capillary  leak.  This  is  the  glass  tube  projecting 
upwards  in  the  figure.  This  leak  was  adjusted  to  have  a  time  constant  of  3 
seconds  (ten  times  longer  than  the  pressure  signal  time  constant).  The  cali¬ 
bration  photo  is  chjwn  in  Figure  7. 

The  two  stopcocks  on  the  side  arms  of  the  sample  chamber  permit 
the  introduction  of  samples. 

The  end  windows  of  the  sample  chamber  and  the  windows  of  the  acous¬ 
tic  chamber  are  set  at  Brewster's  angle  to  minimize  reflection  loss  and  to  im¬ 
prove  the  accuracy  of  calorimeter  readings  of  the  laser  pulse  energy. 

experimental  procedure 

The  sample  was  prepared  by  pssslng  sir  saturated  with  water  vapor 
at  room  temperature  through  the  .ample  tube  for  approximately  15  minutes.  Then, 
the  sample  tube  wa.  clo.ed  except  for  the  capillary  leak.  A  few  drop,  of  water 
remained  In  the  sample  tube  (not  In  the  path  of  the  bam.)  to  lnaure  that  the  air 
would  remain  saturated  with  water  vapor.  The  partial  water  vapor  pressure  was 
approximately  17mm  Hg. 
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Pressure  Step:  -10mm  Hg 
Time  Scale:  1  sec/cm 


Figure  7.  Pressure  Step  Response  of  Sample  Chamber 
With  Capillary  Leak  Open 
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The  ruby  laser  was  pulsed  at  1-minute  20-second  intervals.  A 
calorimeter  measured  the  pulse  energy.  An  oscilloscope  attached  to  the 
Baratron  output  recorded  the  pressure  signals.  The  etalon  and  camera  system 
recorded  the  laser  pulse  wavelength.  A  Fenwall  temperature  controller  and 
thermistor  were  used  to  measure  the  ruby  temperature.  This  temperature  was 
recorded  periodically.  The  coolant  temperature  rose  on  the  average  of  0.2 
to  0.3  of  a  degree  centigrade  for  each  laser  pulse.  After  54  pulses,  the 
temperature  had  risen  from  16.3°C  to  37°C.  The  change  in  ruby  temperature 
caused  the  wavelength  to  scan  over  a  1-1/2  wave  number  range. 

A  test  was  made  to  determine  any  spurious  effects  on  the  signal 
due  to  acoustic  noise  generated  by  pulsing  the  laser.  The  pressure  signal 
was  recorded  when  the  laser  was  pulsed  but  a  stop  prevented  the  laser  beam 
from  entering  the  sample  chamber.  In  each  test  case  no  pressure  signal  was 
observed  on  the  Baratron  output.  Therefore,  the  absorption  measurement  is 
not  disturbed  by  the  noise  of  the  pulsed  laser. 

Timing  signals  were  provided  to  start  the  sweep  in  the  recording 
oscilloscope  and  then  to  trigger  the  laser  pulse  5  milliseconds  later.  This 
system  permitted  observation  of  the  pressure  in  the  sample  chamber  just  be¬ 
fore  the  laser  pulse  and  for  approximately  120  milliseconds  after  the  laser 
pulse.  The  laser  pulse  itself  lasted  for  about  1  millisecond  and  started 
approximately  1  millisecond  after  it  was  triggered.  Thus  by  the  time  the 
first  7  milliseconds  elapsed,  the  laser  pulse  was  over. 

Typical  data  are  shown  in  Figure  8.  The  Baratron  and  oscilloscope 
scale  factors  were  set  so  that  a  1-centimeter  pressure  signal  corresponds  to  a 
1.5H  Hg  observed  pressure  rise.  The  sweep  rate  was  20  milliseconds  per  centi¬ 
meter.  Theie  photographs  are  reproduced  in  one  to  one  ratio.  The  signals  in 
the  upper  photograph  correspond  to  data  plotted  in  the  left  portion  of  Figure  4. 
The  lower  photograph  corresponds  to  data  plotted  in  the  right  portion  of  Fig¬ 
ure  4.  These  particular  data  are  the  signals  resulting  from  pulse  numbers  6 
through  10  and  41  through  45,  respectively.  The  laser  pulse  energy  and  the 
laser  tem>erature  are  also  indicated.  Note  that  as  the  temperature  became 
greater,  the  pul^e  energy  went  down.  The  photographs  from  the  series  etalon 


13 


Report  No.  8670 


PERKIN  ELMER 


Laser 
Pulse 
Energy 
m  joules 


Pressure  Pulses 
20  msec/cm  x  1.5pHg/cm 


Ruby 

Temperature 

c 


Series  Etalon  Photographs 
Etalon  Spacing  5mm 


340  10 


tl 

»  II 

II 

*  II 

u 

*  IJ 

It 

*  11 

11 

♦  II 

266  45 


242  44 


258  43 


254  42 


260  41 


Pulse 

No. 


Figure  8.  Typical  Data  for  the  H_0  Absorption  Spectrum 
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system  are  Bide  by  aide  with  the  photographs  o £  the  pressure  pulses,  end  cor¬ 
respond  to  one  enother  as  lebeled.  the  spacing  of  the  etslon  which  produced 
the  large  outer  rings  was  5  mo. 

DATA  REDUCTION  AND  CALCULATIONS 

A  table  giving  the  laser  pulse  number,  the  laser  puise  energy  reading, 
the  pressure  signal  amplitude,  the  wavelength  shift  data  measured  from  the  series 
etalon  photographs,  and  the  temperature  was  prepared.  Pressure  signal  amplitude 
was  measured  from  the  base  line  established  by  the  start  of  the  pressure  signal 
to  the  peak  amplitude  at  about  1-1/2  centimeters  from  the  start  of  the  trace. 

The  reading  from  the  series  etalon  photographs  was  a  measurement  of  the  diam¬ 
eter  of  the  inner  pair  of  rings  from  the  5  mm  etalon.  When  this  reading  was 
made,  an  average  of  the  intensity  in  each  separate  mode  was  made.  The  diameter 
recorded  is  the  diameter  corresponding  to  the  maximum  intensity  in  the  pulse. 

The  rings  from  the  long  etalon  appearing  as  a  small  spot  in  the  center  of  the 
series  etalon  photographs  were  not  used  since  the  laser  was  not  Q-switched  and 

many  modes  were  mixed. 

The  true  pressure  rise  is  defined  as  the  measured  pressure  rise 
divided  by  a  volume  reduction  factor.  The  laser  pulse  illuminates  only  a 
small  volume  of  the  gas  contained  within  the  sample  chamber.  However,  the 
pressure  rise  is  not  measured  until  it  is  transmitted  uniformly  throughout 
the  volume  of  the  gas.  The  volume  reduction  factor  is  the  ratio  of  the  illumi¬ 
nated  gas  volume  to  the  total  volume  of  gas  contained  in  the  tube.  The  gas 
absorptivity  is  now  calculated  from  the  following  formula: 

a  =  FAP  «r2  pCvT/EP 

where  a  is  the  absorption  coefficient,  AP  i»  the  measured  pressure  rise,  P 
is  the  volume  reduction  factor,  r  is  the  beam  radius,  p  is  the  gas  density, 

C  is  the  heat  capacity  of  the  gas  at  constant  voluma,  T  is  the  absolute 
temperature  of  the  gas,  E  is  tha  pulse  energy  and  P  is  the  nominal  gas  pressure. 
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PLANS  FOR  THE  FUTURE 

The  present  C02  laser-excited  spectrophone  is  somewhat  limited  in 
sensitivity  by  the  long-term  ramp  response  noted  earlier.  Hie  source  of  this 
spurious  signal  must  be  found  and  reduced.  Synchronous  chopping  and  detection 
will  be  used  to  increase  the  sensitivity. 

A  new  sample  cell  for  the  C02  laser  spectrophone  will  be  constructed. 
The  cell  will  be  evacuable  and  capable  of  maintaining  a  pressure  as  low  as  lOp 
Hg.  This  will  make  it  possible  to  measure  gas  absorptivities  at  low  pressures. 

If  it  is  feasible,  we  will  obtain  some  profile  data  on  the  absorption  bands  by 
scanning  the  wavelength  of  the  C02  laser  over  one  of  its  operating  lines. 

We  Will  determine  maximum  sensitivity  for  two  different  cases.  The 
first  case  will  be  sensitivity  to  strong  absorbers  at  low  pressures.  The 
second  case  will  be  the  measurement  of  the  absorptivity  of  pure  helium  contami¬ 
nated  with  known  amounts  of  an  absorbing  gas. 

An  attempt  will  be  made  to  measure  the  modulation  of  the  Stark 
shift  in  anmonia  or  methyl-ethylamine  or  other  molecules  capable  of  strong 
Stark  shift. 

If  time  and  funds  permit,  we  will  observe  the  band  profile  of  water 

O 

vapor  under  self-pressure  broadening  for  the  three  absorption  lines  at  6942. 15A, 
6942. 37A,  and  at  6943. 80^. 


\ 
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APPENDIX 

THE  ABSORPTION  CELL  PRESSURE  RESPONSE 
TO  A  GAUSSIAN  DISTRIBUTED  SOURCE  BEAM 

The  average  pressure  rise  in  the  absorption  cell  can  be  calculated 
from  the  temperature  distribution  and  the  ideal  gas  law.  In  this  analysis  the 
gas  absorptivity  may  be  determined  from  the  average  pressure  rise  if  the  ther¬ 
mal  conductivity  is  known.  The  thermal  diffusivity  can  be  found  from  the 
initial  rate  of  pressure  rise. 

COORDINATES  AND  GEOMETRY  OF  THE  HEAT  FLOW  PROBLEM 

Assume  the  cell  is  long  and  thin,  so  heat  loss  through  the  end 
windows  is  small  compared  with  heat  loss  to  the  walls.  Let  the  laser  beam 
have  a  Gaussian  intensity  distribution  and  cy’indrical  symmetry.  This  re¬ 
duces  the  problem  to  one-dimensional  heat  flow. 

Lot  a  a  coll  radius 
L'(p)  -  beam  intensity  distribution 
I'(p)  =  (W/jf  w2)  exp  <-p2/w2) 
where  W  =  beam  power, 

w  «  beam  spread  (half-width  at  e_1  points). 

Let  r  =  p/a  =  dimensionless  radial  coordinate. 

Thus  the  dimensionless  cell  radius  is  1. 

Let  b  »  w/a,  the  fraction  of  the  cell  radius 
occupied  by  the  beam. 

Thus  I'(r)  =  (W/«  a2b2)  exp  <-r2/b2) . 
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The  intensity  at  the  cell  walls  is  (W/«a2b2)  exp  (-1/b2).  It  is 
desirable  to  have  an  intensity  distribution  which  becomes  0  at  r  *  1.  Hence, 
define  a  nr.vr  intensity  distribution. 

I  (r)  «  (W/ita2b2)  [exp  (-r2/b2)  -exp  (-1/b2)"j. 

When  b  «1,  the  intensity  is  only  slightly  reduced,  but  the  heat  flow  problem 
is  made  much  more  tractable  by  eliminating  infinite  gradients  at  the  tube 
walls. 

Let  a  «*  gas  absorptivity  per  unit  length.  Assume  a  <  <  1,  so  the 

absorbed  power  per  unit  length  is  Cl  W . 

Let  A  =  rate  of  heat  generation  per  unit  volume  per  unit  length 

2 

in  the  cell.  Let  the  unit  volume  be  a  *1  unit  length,  9o 
A(r)  =  (oW/rtb2)[exP  <-r2/b2)  -exp  (-1/b2)]. 

SEPARATION  OF  THE  HEAT  FLOW  PROBLEM 

Let  0  -  0(r, t)  ■  temperature  in  the  cell  above  ambient  tempera¬ 
ture. 

The  heat  flow  equation  is 
b2  9/br2  +  r_1  8  0/ar  -  aVWat-A(r)/K 
where  k  =>  gas  thermal  diffusivity  -  K/pCy 
K  =  gas  thermal  ccoductivity 
p  s  gas  density 

C  =  gas  specific  heat  capacity  at  constant  volume 
v 

Initial  condition:  t(r,o)  ■  0 

Boundary  condition!;  0(1, t)  «  0,  0(o, t)  /  • 
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The  first  boundary  condition  is  justified  because  the  cell  walls 
have  a  much  higher  heat  capacity  than  the  gas.  The  second  boundary  condition 
is  required  on  physical  grounds. 

The  heat  flow  problem  may  be  separated  into  the  steady-state  solu¬ 
tion  and  the  transient  solution.  Let  8  =  u(r, t)  +  v(r) .  Then  v  satisfies 
the  steady-state  equation 

d2v/dr2  +  r_1dv/dr  o  -A(r)/K 

subject  to  v  (o)  t  v(l)  =  0.  Also  u  satisfies  the  transient  solution  equa¬ 
tion 

d2u/dr2  +  r”*"du/dr  =  a2*  ^du/dt 
subject  to  u  (o,  t)  J-  mf  u(l,  t)  =  0 

and  u  (r,  o)  =  -v(r)  .  Obviously  the  steady-state  problem  must  be 

solved  first. 

THE  STEADY-STATE  TEMPERATURE  DISTRIBUTION 

Multiplying  the  steady-state  equation  by  r,  we  have 

d  (r  d  v/dr)/dr  *  -(oW/Krtb':)t^exp(-r2/b2  -exp  (-l/b2)J. 

-  Integrating  and  dividing  by  r,  ve  obtain 

dv/dr  m  (OW/Kx)  £exp(-r2/b2)/2  r  +.exp  (-l/b2)r/2b  J  +  C^/r. 

The  first  exponential  factor  can  be  Integrated  term  by  term.  Replace  it  with 
l/2r  +  (1/2)  )f|l(-l)n  r2nml/b2nnl 
and  integrate  to  obtain 

v  -  (oW/itK)[(l/4)  %ml  (-l)n(r/b)2n/nn! 

+  (1/2)  An  r  +  (r2/4b2)  exp  (-1/b2)]  +  Cj  in  r  +  Cj. 

To  satisfy  v(o)  t  •,  let  «  -aW/2kx 
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To  satisfy  v(l)  a  0,  choose  so 
V  -  (»M«K)[C1(-l)”('2”-l)/»2”nn!  *  exp  (-1/b2) (r2-l)/b2] 


THE  STEADY-STATE  PRESSURE 


Tho  relation  between  the  pressure  rise  and  the  temperature  rise  is 
found  by  differentiating  the  ideal  gas  law, 

-  dT/T  where  P  is  the  nominal  pressure  and  T  is  the  ab¬ 
solute  ambient  temperature.  The  average  pressure  p  is  found  by  the  integral 

P1  =  (2it  1  unit  length/V)  (P/T)  f1  v  (r)  r  dr 

*c 

where  the  dimensionless  volume  V  -  n  l2  1  unit  length.  Hence 

PX  =  (2P/T)  f1  v(r)  r  dr 
Jo 

=  (QWP/4itKT)  {Cl<-»1>2n+2/(„+l)  -  r2]  /b2nn„! 

+  exp  (-l/b2)(r4/2-r2)/b2\1 

Jo 

"  <aWA«r)[Ci<-1)n+l/b2n(»fl)!  -  exP(-l/b2)/2b2] 

"  <aWP/4*KT)[l-b2  +  exp(-l/b2)(b2-l/2b2)] 

Let  B«  B(b)  a  [l-b2  +  exp(-l/b2)(b2-l/2b2)] 

Then  the  absorptivity  is 


a  «  P24ttKT/WPB 


(A-l) 


Table  A-l  gives  values  of  B  for  cases  of  interest. 
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TABLE  A-l 

b  a*  0  0.20  0.25 

B  1  0.960  0.938 

THE  TRANSIENT  TEMPERATURE  SOLUTION 

A  solution  of  the  transient  temperature  solution  yields  the  gas 
thermal  diffusivity  in  terras  of  the  amplitude  and  initial  slope  of  the  pressure 
rise . 

Assume  that  u  (r,t)  =  R(r)  T(t).  The  heat  flow  equation  becomes 
R_1(d2R/dr2  +  r"1  d  R/dr)  •  a2*-1!-1  dT/dt  -  -?2 

2 

where  -?  is  an  arbitrary  constant. 

The  solutions  are 

2  2 

T  «  exp(-K  5  t/a  )  and 

R  .  C  J  (?r)  +  DN  (?r) 

O  o 

Since  u(o,  t)  +  *,  D  =  0.  To  satisfy  u(l,t)  a  0,  chose  ?  a  where 

is  the  mth  root  of  J  (?)  =  0. 
hi  o 

The  general  solution  so  far  is 

u<r.t>  -  Cl  Jo  «V>  ^  t/,2) 

To  satisy  the  initial  condition, 

u(r,o)  -  Cl  V.  "«r)  -  'T(r)’ 

nulciply  both  sides  of  the  equation  by  rJ0<?kr),  integrate  iron  0  to  1,  and  use 
the  orthogonality  properties  of  the  Bessel  functions  to  obtain 


0.30 

0.910 
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Ck  -  -r 2/J12<^)](o«/4JtK)  - 

•[rn-1(-l)n(n'.nb2n)"1  J  (r2n+1-r)Jo(?kr)dr 

1  0 

+  b  2  exp(-l/b2)  J  (r3-r)Jo(«;kr)drJ 
-  [2/J^Sjj)  Sk  ]  (QWMnK)- 

•{Ci<-i>n<n!nb2n)’1  [l-  Cvji<  v)  f  >+1w-] 

+  4  exp  <-l/b2)/b2?k2  } 

THE  TRANSIENT  PRESSURE  SOLUTION 

The<  transient  pressure  may  now  be  determined  as  before. 

P2(t)  -  (2P/T)  Jq  u  (r)  r  dr 

«  (2P/T)  E*  C  J.  (?)  •~l  exp  (-H  ?  2  t/a2) 

m«i  m  1  m  m  r  m 

=  (QWP/4jacn  E*  .  D  exp(-H  ?2  t/a2) 

Dbil  IQ  IQ 

where  D  »  (4/?2  )  U  exp(-l/b2)/b2?2 

tq  m  m 

*  0-»”  [»-  (yw)  f  t2ntlj„<v>dr}'':"b2''} 

o 

The  Integral  in  Dm  decreases  with  Increasing  n  and  with  Increasing 
m.  It  may  be  evaluated  by  successive  applications  of  integration  by  parts. 

We  found  evaluation  by  direct  Simpson  integration  to  be  simple  and  precise. 

Table  A-2  gives  the  coefficients  for  ss.1  through  20,  and  for  various  values  of  b. 
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TABLE  A- 2 

VALUES  OF  D 

m 

b 

0 

0.20 

0.25 

0.30 

m 

1 

-1.1080 

-1.0358 

-1.0123 

-0.9728 

2 

0.1398 

0.10A5 

0.8683 

0.070A 

3 

-0.0A55 

-0.0211 

-0.01A1 

-0.008A 

A 

0.0210 

0.0056 

0.002A 

0.0009 

5 

-0.0116 

-0.0010 

-O.OOOA 

-0.0001 

6 

0.0072 

0.000 A 

7 

-0.00A8 

8 

0.003 A 

9 

-0.0025 

10 

0.0019 

11 

-0.0015 

12 

0.0012 

13 

-0.0010 

1A 

0.0008 

15 

-0.0007 

16 

0.0005 

17 

-0.0005 

18 

0.000A 

19 

-0.000A 

20 

-0.0003 
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THE  TOTAL  PRESSURE  RISE 

Let  p  =  total  pressure  rise  =  p^  +  p^t), 

P  =  <0WP/4TtKT)  [B<b)  +  if  .D  exp  (-H  ?2  t/a2)l 
l  m=l  m  r  m  'J 

Let  P(»)  =  A,  the  final  amplitude  of  the  step  response. 

The  diffusivity  can  be  obtained  from  the  amplitude  and  the  initial 
slope.  Initial  slope  =  p(o) 

=  -( JWP/4 JtKTjKa"2^00  ,D  P2 
m=l  m  m 

=  -  iWP/4jrKT)  B(b) _ 


A _ 

“  T  C(b) 

where  t  is  the  principal  time  constant 


2  2 

a  / (2.4048)  h,  and  C(b)  is  a  correction 
time  constant.  The  correction  factor 

factor 

to  give  the  apparent 

is  given 

C(b)  =  -B(b)/E*  D  (?  /?  )2 
m=l  m  m  1'  ' 

in  Table  A-3  for  various  values  of  b. 

TABLE  A-3 

b 

~0  0.20 

0.25 

0.30 

C 

1.48  1.44 

1.35 

1.32 
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CONCLUSIONS 


where 


where 


The  absorptivity  is  given  by 

a  =  4rtKTA/WPB,  (A-2) 

K  =  gas  thermal  conductivity 
T  =  ambient  gas  temperature 
A  =  final  amplitude  of  step  response 
W  =  beam  power 
F  =  nominal  gas  pressure 

B  =  a  function  of  b,  the  ratio  of  beam  radius  to  cel}l  radius, 
given  in  Table  A-l. 


The  thermal  diffusivity  is  given  by 

K  =  p  (o)  a2C/ (2 .4048) 2  A, 

p(o)  =  the  initial  rate  of  pressure  rise, 
a  =  cell  radius 

C  =  a  function  of  b  gi-'en  in  Table  A-3. 


(A-3) 
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